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Abstract 



We extend the now classic structural credit modeling approach of Black and Cox to a class 

of "two-factor" models that unify equity securities such as options written on the stock price, 

^ and credit products like bonds and credit default swaps. In our approach, the two sides of the 

^ stylized balance sheet of a firm, namely the asset value and debt value, are assumed to follow 

or* a two dimensional Markov process. Amongst models of this type we find examples that lead 

to derivative pricing formulas that are capable of reproducing the main features of well known 

equity models such as the variance gamma model, and at the same time reproducing the stylized 

facts about default stemming from structural models of credit risk. Moreover, in contrast to one- 

factor structural models, these models allow for much more flexible dependence between equity 

OO and credit markets. Two main technical obstacles to efficient implementation of these pricing 

formulas are overcome in our paper. The first obstacle stems from the barrier condition implied 

Q by the non-default of the firm, and is overcome by the idea of time-changing Brownian motion 

^ in a way that preserves the reflection principle for Brownian motion. The second obstacle is the 

y—i difficulty of computing spread options: this is overcome by using results in recent papers that 

J> make efficient use of the two dimensional Fast Fourier Transform. 

• ^ 

X 

Key words: Credit risk, capital structure, structural model, Levy process, time changed models, 
equity derivative, implied volatility, credit spread, spread options. 



1 Introduction 

Merton [20|, Black and Cox and other pioneering researchers in credit risk well understood 
that dynamics of a firm's equity and debt should be modeled jointly and that credit derivatives and 
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equity derivatives are linked inextricably. To this day, however, it has proved difficult to capture the 
dynamical essence of these two aspects of a firm's capitalization. The papers by Leland [16] and 
Leland and Toft [fTTll provide a conceptual basis, but they remain strongly attached to the framework 
of diffusion processes and have a one dimensional source of randomness . 

The above structural models can all be classified as one-factor models with the asset process as 
the only underlying source of randomness. Such models have the severe limitation that the firm's 
equity and debt are perfectly correlated (i.e. they are related by a deterministic function), while it is 
clear in real life that firms, to a greater or lesser extent, have stochastic liabilities that are not perfectly 
correlated with assets. As an extreme illustration, hedge funds, with their long/short positions, typi- 
cally have liabilities as volatile as their assets. Thus the fact is clear that to accurately model capital 
structure, a stochastic factor process of dimension at least two is necessary. In the context of contin- 
uous time finance, the technical and computational challenges implied by this fact have not yet been 
satisfactorily addressed, and these challenges are the main focus of the present paper. Only a few 
authors have been able to make substantial headway in modeling actual observed capital structures 
by two factor models. Eberlein and Madan [|TOl , in a recent working paper, have treated firm asset 
and liabilities as imperfectly correlated processes. Equity equals the asset/liability spread, and they 
use this fact to calibrate both the asset and liability values from the observations of a firm's implied 
equity volatility surface. 

A second deficiency in the standard structural framework is the reliance on diffusion processes, 
with the consequence that default events are predictable and so instantaneous default is either certain 
or impossible [21 J. Thus in such models the short spreads are either infinity or zero, counter to 
the fact that short spreads are observed to be positive even for investment grade firms. The natural 
way to overcome this deficiency is to introduce jumps into the asset process. A number of authors, 
notably [|ll,[l4l, have successfully implemented jump diffusion and pure jump versions of the Merton 
model. However they share that model's unrealistically simple debt structure. Similar extensions to 
the Black-Cox first passage framework, however, have had only a limited success, due to the technical 
difficulty of solving the first passage problem for jump processes. The Kou-Wang model [15J with 
exponentially distributed jumps was able to work because of the special nature of the underlying 
process. 

This difficulty with predictable defaults was the original motivation for replacing structural models 
by reduced form models [11411 and incomplete information models [191. Recently, a class of "hybrid" 
reduced form models that include the stock price and a default hazard process have been developed. 
These model equity and debt products more realistically by allowing the stock price to jump to zero at 
the time of default. Carr and Wu [7J take a stochastic volatility model for the stock price and assume 
that the default arrival rate is driven by the volatility and another independent credit risk factor. In 
Carr and Linetsky [|6l, the stock price has a local volatility with constant elasticity of variance, and 
the default intensity is specified as an affine function of the instantaneous variance of the stock. [19J 
obtain even more generality while retaining analytic tractability by applying a time change to the 
hazard rate and stock price processes. All three frameworks are able to capture the so called leverage 
effect and the co-movement between volatility and default intensity. 

Within the structural credit modeling framework, [flTl generalizes the Black-Cox model by treat- 
ing the log-leverage ratio Xt := logVt/ K{t) as a time-changed Brownian motion (TCBM), where 
e~^^Vt := e^* denotes the (per-share) firm asset value process discounted at a constant rate r and K{t) 
is a deterministic default threshold. The time of default is the first passage time for the log-leverage ra- 
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tio to cross zero. Like other structural approaches along these lines, this model cure the inconsistency 
with observed short spreads and add the flexibility to include jumps to default and volatility clustering. 
One contribution of the TCBM approach in [1 IJ lies in an innovative mathematical treatment of the 
first passage to default that allows the reflection principle and corresponding first passage formulas 
for Brownian motion to extend to a broader class of processes, leading to analytical tractability in a 
more general setting. 

The object of the present paper is demonstrate how to embed the Black-Cox framework in a simple 
way into a two-factor framework that allows the firm's equity and debt to be partially correlated while 
retaining tractability of the underlying default model. We do this by treating the default threshold K 
as a positive stochastic process Dt := e^*"*""^* that we can think of as the market value of the firm's 
liabilities, per share. Put another way, we treat the firm's debt or liability as a new stochastic process, 
not fully correlated with the asset process Vt. If we consider as well the stock price St (assuming the 
number of shares is constant and the firm pays no dividends) and the log-leverage Xt, then a minimal 
set of additional assumptions for combining these processes is: 

Assumption 1. The pre-default dynamics of any two of the four processes Vt, Dt, St, Xt is Markovian 
and time-homogeneous, and determines the dynamics of the remaining two processes by the equations 

St = max{Vt-Dt,0), Xt = \ogVt/Dt. (1) 

We assume that the discounted processes e^^^Vt, e~^^Dt, e~^^St are martingales under some risk neu- 
tral measure Q, and the interest rate r is constanj^ The time of default is 

f = mi{t\Xt < 0} = inf {t\St = 0}, (2) 

and after default Xt = St = 0. At the time of default, all securities are assumed to be valued in terms 
of a "recovery" random variable R. 

In this paper, we will make additional restrictive assumptions on the form of V, D,S,X to obtain 
a workable tractable framework. Under these restrictions, the pure credit dynamics of two factor 
models with a constant recovery will be consistent with the TCBM credit framework of [1 IJ. In [ 13] 
good calibrations of such credit models to a large dataset of CDS spreads for Ford Motor Company 
were obtained, thus verifying the quality of the framework as a model for credit risk. In addition our 
two-factor models price equity options as barrier spread options on V, D. Thus option pricing in two 
factor models faces the type of computational challenges for spread options that have been studied in 
such papers as [8^, f5\ and We are able to use the TCBM structure and properties of Brownian 
motion to develop an efficient equity option pricing algorithm that uses a two dimensional Fast Fourier 
Transform (FFT). Such a fast algorithm is needed not just in the forward direction for security pricing, 
but more importantly to solve the "inverse problem" that arises in calibration to a dataset of observed 
security prices. In this paper, because we have efficient pricing of the basic securities to be used in the 
calibration, we are able to demonstrate the feasibility of efficient statistical estimation of two factor 
models to a dataset of CDS curves and implied equity volatility surfaces. 

The above assumptions on the firm's capital structure can only be valid for a period [0, T] over 
which the firm pays no dividends or debt coupons, and does not issue new shares or debt. A consistent 
firm model that incorporates such real world features will also be of interest for future research. 

'While the arbitrage pricing theory requires only that e~''*5f be a martingale, we make a stronger assumption to 
simplify the framework. 
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This paper will discuss several implementations of the two-factor framework where 5*, D and V 
are discounted martingales. In all these implementations, we are able to overcome the two technical 
obstacles, namely the treatment of the first passage to default and the efficient computation of spread 
options. In §2 and §3 we investigate the case where Vt and Dt are correlated geometric Brownian 
motions (GBMs). The resultant default model extends the Black-Cox model and shares its known 
shortcomings, such as zero short spreads. As well it tends to generate rather unrealistic implied 
volatility surfaces. Therefore in §4 we allow Vt and Dt to be Levy subordinated Brownian motions 
(LSBMs) driven by a single time change, in this case either a gamma process or a process with 
exponential distributed jumps. We investigate some of the possible shapes of both the CDS curve 
and the implied vol surface. In §5, we investigate how to calibrate such models to market CDS and 
implied vol data on a single date. We then exhibit the results of a simple calibration of the GBM 
and LSBM models to data for a typical default risky firm. Ford Motor Company. Finally, §6 offers a 
summary and some directions for future exploration. 

2 Risk-Neutral Security Valuation 

As we have explained, the firm's capitalization is modeled by the four processes V^i, Dt, X^, St satis- 
fying Assumption ([T]), and default t* is the first time Vt < Dt, or equivalently when St = 0. We work 
in a risk-neutral filtered probability space J-', {J-'t)t>o, Q), where Vt, Dt, Xt, St are adapted to the 
filtration J't and t* is an J't stopping time. If Qt C J-'t denotes the "market filtration", note that the 
stock price St is Qt measurable, whereas X, V and D are not. In practice Xt,Vt, Dt must be inferred 
from the prices of securities trading on the firm, and possibly its quarterly balance sheets. 

In this section, we demonstrate that the price of basic credit and equity derivatives can be reduced 
to computations involving the joint characteristic function of Yt := [vt, dx] conditioned on non- 
default at time T: 

^Mui,U2;T;vo,do) := <,,Je^^"^''-+^^^-)l|t.>T}] • (3) 
As a special case, note that the probability of no default at time T is: 

P{T;vo,do) = ^ND{0,0;T;vo,do) = E^, ,Jl{t*>T}] . (4) 

The reader can anticipate that in subsequent sections, we will introduce a number of models of 
dynamics where a generalized reflection principle holds and implies that computations involving ^nd 
can be reduced to computations involving the unconstrained characteristic function 

^u,,U2;T;vo,do) := E^^^,Je^(-^^-+"^'^-)] . (5) 

Moreover, our models will have a common property on Yt. 

Assumption 2. For any t > 0, the increment Yt — Yq is independent of Yq. 

This implies that the characteristic function $ of Yt factorize^ 

Ero[e™^^] = e™^o'$(w; T), T) := Eyje*"^^^-^")'] . (6) 

where T) is independent of Yq. Thus, in the more specialized setting we have in mind, pricing 
of all important derivatives will be reduced to computation of explicit low dimensional integrals. 

^Here and subsequently we adopt matrix notation M = [7/1,^2], = [ut, c?(], and in particular F' denotes the transpose 

of y. 
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2.1 Defaultable Bonds and Credit Default Swaps 



At any time t prior to default t*, we consider the value of a zero coupon bond that pays $1 at maturity 
T if the firm is solvent. In the event that t* < T, it might be reasonable to suppose that the recovery 
value of the bond will be dependent on Dt*. However, to avoid a detailed analysis of stochastic 
recovery modeling, we make a mathematically simple hypothesis: 

Assumption 3. The recovery value of a zero coupon bond with maturity T, at the time of default 
t* < T is a constant R e [0, 1). 

This assumption is analogous to the recovery of par mechanism often made in credit risk mod- 
eling, and it shares some of its limitations. Then one has the risk- neutral valuation formula for the 
pre-default price Bt{T) at time t of the maturity T zero coupon bond: 

Bt{T) = [e-^(^-*)l|,,>T} + i?l|,<,.<T}e-^(**-*)| J-,] , (7) 

which leads to: 

Proposition 1. 1. The pre-default price at time t of a zero coupon bond with maturity T with 
recovery of par is given by 

Bt{T) = e-^<^-*)p(r - t; Vt, dt) + Re-''^^~'\l - P{T - t; Vt, dt)) . (8) 

2. The fair swap rat^\for a CDS contract with maturity T = NAt, with premiums paid in arrears 
on dates tk = kAt, k = 1, . . . , N, and the default payment of (1 — R) paid at the end of the 
period when default occurs, is given by 

(1 - R) fcf=i'[l - Pih; vo, do)] [e-"-'" - e-*^+^] + e-'^[l - P(T; Vo. d^)] 

CDSiT] Vo,do) = Tf 

AtZtiPitk;vo,do)e~^^^ 

(9) 



2.2 Equity derivatives 

We have assumed that St = for all t > t*. This is a plausible idealization of the observed fact that 
stocks typically trade near zero for a period after a default plus the fact of limited liability that ensures 
S't > 0. By the martingale assumption it follows that for any t < s prior to default 

S, = E\-'^'-'\Vs - Ds)l{t,>s}\J't] = (Vt - A)l{t.>t}. (10) 
The second equality comes from Doob's optional stopping theorem. We notice that 



sAt* 



is a Q martingale evaluated at a bounded stopping time s At*, which is also a Q martingale. In ( [TO] ), 
St is independent of the debt maturity s. This is different from the standard setup in the Merton model 

^CDS contracts were restructured at the "big bang" on April 8, 2009, and since then quoted CDS spreads must be 
reinterpreted. We ignore this complication here, but see |2| for a discussion. 
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and Black-Cox model, which makes it more parsimonious. Moreover, the time t price of a maturity 
T > t forward contract with strike K will be St — e^'''^^^^^K. A European call option with (positive) 
strike K and maturity T has time t pre-default value 

Callf^ = E^[e-^'(^-*)(rT -Dt- KYl{t*>T}\J't]. dD 

Observe that this is equivalent to a down-and-out barrier spread option with a leverage barrier on 
Xt = 0. Put-call parity also holds in such a model, implying that Califs — Putf ^ = St — Ke^^^'^^^K 
When a closed or computable form exists for the non-default characteristic function ^nd the 
above option pricing formula is amenable to Fourier analysis, following the method developed in [[T2l| 
for vanilla spread options. There it is proved that the spread option payoff function has an explicit 
two-dimensional Fourier transform: 

Proposition 2. For any real numbers e = (ei, £2) with €2 > and ei + €2 < — 1 

- - 1)+ = (27r)-2 / / e'("i"^+"^"2)p(ui,M2)ci'?x (12) 



=^1 



'+it 



where P{ui,U2) = ^(^("i+^^MM ^U2) ^ where T{z) is the complex gamma function defined for > 
by the integral 

poo 

Viz) = / e-H'-^dt . 



CalF(t;o, (io) = 77^^ / / ^NDiui,U2;T;vo, do)P{ui,U2)d^u . (13) 



Combining this formula with the Fubini Theorem leads to the following formula for a call option 
with strike K = 1 and maturity T: 

-rT 

{2V) 

For a general strike K = e'^, we use homogeneity to write 

Call^^(t;o, do) = KCalf{vo -k,do-k). 

Such explicit double integrals are sometimes efficiently computable for a full range of vq, do values 
using a single two-dimensional Fast Fourier Transform. 



3 Geometric Brownian Motion Hybrid Model 

We now consider the two factor model where Vt = e^'*"^^' , Dt = e''*+'^' are jointly given by a two- 
dimensional geometric Brownian motion: 

^ =rdt + a^dWt, — = rdt + aadZf, dWtdZt = pdt. (14) 
Vt Dt 

In this case, the stock price St = Vt — Dt and log-leverage ratio Xt = vt — dt follow the SDEs 

dSt dVt - dDt 



St Vt - Dt 

1 

2' 

Intuitively one normally expects to find ciy > ad> 0, and to have the correlation p E (—1, 1). 



(15) 



dXt = --{al - a1)dt + a^dWt - addZt . 
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3.1 Stochastic Volatility Model 

Before investigating the form of $ and hence the pricing formulas, it is worthwhile to note that this 
two factor model is identical to a specific stochastic volatility equity model, analogous to the Heston 
model. To see this, first we note that we can write 

Xt^Xo + ax[at + Bt] (16) 

2 2 

Here — cr^ — 2 pay a a + cr^ and a — '^^ , and the Brownian motion B is correlated to W, Z with 

dBdW = Pvxdt, GxPvX ^CTv- craP 
dBdZ = Pdxdt, axPdx = pcTv - crd ■ 

Next we apply the Ito formula to obtain the SDE for the pre-default stock price 

dSt _ r(e^* - e'^*)dt + {a^e'^^dWt - aae'^'dZt) 

St ~ e'"^ - 

(Tye^'dWt- aadZt 

= rdt^ . 

e^* - 1 

The martingale term has stochastic quadratic variation with increment dS^/S^ — f{Xt)dt where 



m feCi' (17) 



Furthermore, the cross variation increment is dXtdSt/St — g{Xt)(Tx-\fJ{Xt)dt where 



ale" - p(Tv(7d{e'' + l) + al 
(Jx^/{(Tve'= - (Tdpy + (1 - p'^Wd 



g{x) := ^"-"-^^^ ^^>^-<i . (18) 



Therefore, using the Levy theorem to define a new Brownian motion, one can prove 

Proposition 3. In the GBM hybrid model, there are independent Brownian motions B,B^ such that 
the log-leverage process is given by 

Xt^Xo + ax[at + Bt] 
and the stock price follows a stochastic volatility process 

dSt/St = rdt + at[psx,tdBt + psx,tdBl] (19) 

with 



at = fiXt), psx,=9{Xt), psx,t = V'^-9iXtr- (20) 
Moreover, the default time t* is the first passage time for Xt to cross zero, and is predictable. 
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Remark 1. The processes Vt, dt can be expressed in terms of the independent drifting BMs Bt + at 



and + aH where a-^ = = 



Vt = Vo + a^[pvx{Bt + at)+p^x{B^ + aH)] (21) 
dt = do + ad [pdx{Bt + at) + pdx{B^ + aH)] 



where 

pvx = a/1 - p^dd/ax, pdx = a/I - p^o-ujox ■ 
Finally, we note that in the GBM hybrid model, the explicit characteristic function is 



$^^^(«;r,ro) = exp ^uY^--uLv^ -i — {alal) 



(22) 



where E = [a^, pOyOd; pa^ad, aj] 



3.2 Pricing 

Basic securities we need to price, namely, defaultable bonds and equity call options, have payoffs 
that vanish on the set {vt < dx} and are subject to a "down-and-out" barrier condition. The next 
proposition shows how the barrier condition can be easily dealt with for such securities. First we note 
that the linear change of variables [Xt, X^^]' = MYt, Yt := [tn, dt]' for the matrix M — [1, —1; 1, m] 

2 

with m = ^^"^^1^2 leads to independence of Xt — vt — dt and X^. This fact allows us to state and 
prove the following important result: 

Proposition 4. Consider an option with maturity T and bounded pay ojf function F{v,d) that vanishes 
on the set {v < d}. Let f{vo, do; T) denote its value at time 0. In the geometric Brownian motion 
model, the down-and-in barrier option with initial state vq > do and terminal payoff F is equivalent 
to a vanilla option with the same payoff, but with linearly transformed initial state [vq, do] and an 
extra factor. Precisely, 

fDi{vo: do; T) = e-2"(^°-'^°)/'^-/(i;o, do; T) (23) 
where [vt, dt]' = R[vt, dt]', R = M-^[-l, 0; 0, 1]M. 

Proof: Note that the matrix i? is a skewed reflection matrix, which hints that this result is essentially 
a consequence of the reflection principle for Brownian motion. By intermediate conditioning, 

E„o,do[F(^;T,dT)l{r<r}]=E^x[Exo[F(M-i[XT,X^]')l{x,>o}l{r<r}|^^]] ■ 
For fixed X^, the reflection principle the inner expectation can be written as an integral 

where G{x, y) := F{M~^[x, y]'). Here we have used a standard result for Brownian motion condi- 
tioned on crossing a barrier. The vanilla option with the same payoff can be written 
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where the inner expectation equals 



Jo cixVT V (TxVT J 



The desired result follows by a direct comparison of these two formulas for the inner expectations. 

□ 

Corollary 5. In the geometric Brownian motion hybrid model 

1. The survival probability by time T can be written 

P\t* > T\vo,dQ] = P[vt > dT\vo,do\ - 6'^''^°/'''' P[vt > drlvoJo] 
and the price of a zero-recovery defaultable zero-coupon bond is given by e~'"^P[t* > T|i;o, c^o]- 

2. The price of an equity call option can be written 

F{vo, do; T) - e-2"^°/'^-F({}o, d^-, T) 

where F{vo, do; T) = e'^'^KyQ^do [{^^"^ ~ ^'^^ — K)^l!^y^^dT}\ ■ vanilla call option price with 
maturity T and strike K = 1 can be computed by the two-dimensional FFT: 

F{vo, do; T) = / / T; vo, do)P{ui, U2)d\ (24) 



4 Levy Subordinated Brownian Motion Hybrid Models 

We have seen that the two-factor GBM model implies that the stock process 5 is a rather specific 
stochastic volatility process with continuous paths. Moreover the log-leverage process X is an arith- 
metic Brownian motion with constant drift, and the resultant Black-Cox credit model is well known 
to be unable to capture the fine effects in observed credit spreads. 

The time-changed Brownian motion (TCBM) credit framework of [fTTll introduces a non-decreasing 
"time-change" process Gt independent of B and replaces the Brownian log-leverage process by its 
time-change Xt = Xo + ox {Pct + OLGi\ to create a much richer range of dynamics, allowing for 
purely discontinuous components as well as "stochastic volatility". The relevant notion of default 
by first-passage of the log-leverage process to zero has been well understood in this setting. A 
non-decreasing Levy process Gt is called a subordinator, and a Levy subordinated Brownian mo- 
tion (LSBM) in general includes purely discontinuous components. Any LSBM Wct + oiGt has the 
independent increment property and is Markov and therefore we can say it is a one-factor process. An 
important consequence of this one-factor property is that it excludes stochastic volatihty effects that 
by definition involve further factors. 
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The same time-change ideas can be applied to our two-factor GBM hybrid model, and will provide 
a dramatic increase in flexibility to match effects observed in market data. To retain the simplicity 
of two-factor models, we focus here on the LSBM case with a single subordinator Gt and the two 
uncorrelated drifting Brownian motions Bq^ + aGt, Bq^ + a-^Gt- We assume the natural filtration Ft 
contains a{Gu, By, B:^ : < u < t,0 < v < Gt}. 

The assumptions underlying the LSBM two-factor hybrid model are: 

Assumption 4. 1. The time-change process Gt is a Levy subordinator with mean E'5[G'(] = t. 

2. The log discounted firm value Vt = rt + \og{Vt) and log discounted firm liability dt = rt + 
log{Dt) are both LSBMs, with the same time change Gt, i.e. 

Vt = vo + cr^ [pyxiBct + aGt) + pvxiB^, + a^Gt)] (25) 
dt = do + ad [pdxiBct + aGt) + PdxiB^^ + a-^Gt)] • 

Here, the parameters are chosen as in section 3. 

3. The log-leverage ratio Xt := logiVt/Dt) = Xq + axlBct + aGt] is also a LSBM, and St = 
Vt - Dt. 

4. The time of default is t*, the first passage time of the second kind for X to cross zero, defined 
by 

t* = M{t\Gt > t} (26) 
where r = mf{t\Bt + at < — Xq/ctx}- All processes are stopped at t*. 

5. The interest rate r and recovery fraction R are assumed constant. 

In the model calibration that follows in Section |5] we will consider two specific forms for the 
subordinator Gt'. 

1. The first type of time change is an exponential (EXP) jump process with constant drift, that is, 
G has characteristics (6, 0, u) where b E (0, 1) and iy{z) = ce~^l°- ja, a > on (0, oo), the Levy 
measure, has support on M+. The Laplace exponent of Gt is 

^^^P(M,t) := -logE[e-"^*] = t 

and by choosing a = ^ the average speed of the time change is normalized to 1 ; 

2. The second type of time change is a variance gamma (VG) process [|T8l , that is, G is a gamma 
process with drift having characteristics (b, 0, v) where b E (0, 1) and ^{z) = ce^^/"/^;, a > 
on (0, oo), the Levy measure, has support on IR+. The Laplace exponent of Gt is 

^^^{u, t):=- log E[e-"^*] = t[bu + c log(l + au)] (28) 

and by choosing a = ^ the average time change speed is normalized to 1; 



bu + 



acu 
1 + au 



(27) 
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The practical consequence of the precise way the time-change is introduced, and in particular the 
associated definition of default t* as a first passage of the second kind, is that all expectations relevant 
for the pricing of securities can be done efficiently by iterated expectations. For example, we have a 
simple formula for the characteristic function of {vt, cLt)'- 

Since given by (|22]) has the nice feature that the T dependence takes an exponential affine form 

which implies that the GBM pricing formula easily extends to TCBM with a Levy subordinator. 

Proposition 6. Consider an option with maturity T and bounded payoff function F{v,d) that pays 
only ift* > T. Let f'^^'^^VQ, do; T) denote its value at time under the GBM hybrid model, and 
j-LSBM^y^^ (io; T) its value under the LSBM model. Then 

f'^'Hvo,do;T)=E'^[f^'\vo,do;GT)] 

Proof: If we write vt = vgj,, dx = dcj. and define r as in ( |26l ), we know that /'^^*^(fo, c?o; T) = 
Ko.doi^i^T, t^r)l{r>T}]- Then in the LSBM model, 

/^^^^(t;o,cio;T) = EQ[Ej^,jF(5G.,rfG.)l{.>G.}|GT]] (29) 
= E«[/«^^(.;o,cio;GT)]. 



□ 



As an important example, we can see that combining the above result with Corollary [5] leads to the 
following formula for the equity call option with maturity T and strike K = 1 in any LSBM model 
where the time-change G has Laplace exponent ip: 



^LSBM 



{vo, do; T) - e-2"K-rfo)/-x^i-5BAf^~^^ j^. ^^^^ 



where the vanilla spread option price is 



■^LSBM, 



-rT 



vo, do; T) 



exp 



iuYo-tpiuJ^u /2-iu{a^^,aj)' /2,T) P{ui,U2)d^u. (31) 



5 Calibration of LSBM models 

The aim of this calibration exercise is to demonstrate that the simple two-factor LSBM hybrid frame- 
work is capable of fitting simultaneous market CDS and implied volatility prices on a firm, in this 
case Ford Motor Company, at any moment in time. We chose Ford as an example of a large, highly 
traded, firm, that has been very near to default in recent years. We do not here attempt a large scale 
survey of how the model performs for a broad range of firms over different periods of time. However, 
we will see encouraging results from our small study, that suggest that acquiring and analyzing such 
a dataset may be worth the considerable expense and effort involved. 
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5.1 Data 



We observed equity and credit market data for Ford Motor Co. obtained from Bloomberg at two 
moments during the post credit-crunch period: once on July 14, 2010 and once on February 16, 201 1 . 
On these dates we noted: 

1. The stock price was $11.81 and $16.05 respectively; 

2. Midquote implied volatilities /Vd ^ for moneyness 

V := {0.4, 0.6, 0.8, 0.9, 0.95, 0.975, 1, 1.025, 1.05, 1.1, 1.2, 1.3, 1.5} and with times to maturity 
r := {37, 65, 156, 191, 555} calendar days on July 14, 2010 T := {30, 58, 93, 121, 212, 338} 
calendar days on February 16, 201 1; 

3. Midquote CDS spreads CD St for tenors f :={1,2,3,4,5, 7, 10} years; 

4. US treasury yields for maturities 

r := {Im, 3m, 6m, ly, 2y, 3y, 5y, 7y, lOy, 20y, 30?/}. 

From 2009 to 2011, Ford Motor Co. steadily recovered from its near default during the 2007/08 
credit crunch and expectations from the financial market correspondingly rose. This improvement 
in the firm's fortunes is manifested in an observed decrease of both its implied volatilities and CDS 
spreads between the two observation dates. 

Remark 2. We found that deep in-the-money (ITM) options are not very liquid and deep out-of- 
the-money (OTM) options are difficult to control numerical errors as their prices are very close to 



zero. For very short time to maturity options, our FFT formulas ( [24] ) and pT] ) are subject to higher 
truncation errors as the integrand does not decay fast enough. Therefore in our calibration, we did 
not use implied volatility data with extreme moneyness D = 0.4, 1.5 and with short time to maturity 
(TTM) T = 30, 37, 58 calendar days. 

5.2 Daily Calibration Method 

To test the GBM, EXP and VG models described above, we performed independent daily calibrations 
to the data observed on the above two dates, making use of the model formulas (|9]) and ( |30l ). We 



mention that we adjusted these formulas to incorporate the yield curve Yt obtained from a boot-strap 
of the observed US Treasury yields, by the replacement of all discount factors e"*^* with corresponding 
factors e~^**. 

It is natural to assume that stock prices are perfectly liquid and hence are perfect observations of 
the process St = e~^*[e^* — e''*]. On the other hand, CDS and equity option markets are much less 
liquid therefore the N = 7 observed CDS spreads and M ~ 50 observed implied volatilities are not 
assumed to match our model prices exactly. Thus at any moment, St is exactly observed, while Xt 
must be filtered from the market data. 

On any given date t, under these assumptions, the risk neutral parameters for both LSBM models 
to be calibrated are = (p, a^, (Jd, b, c, R, Xt) G M x M^. The GBM hybrid model nests inside both 
LSBM models for any value of b as the limit with c = oo. 

Our method is a relative least-squares minimization of the error between the model and market 
CDS spreads and implied volatilities observed at a single instant in time. For each T E T, let CDSt 
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and CDSt{Q) denote market and model CDS spread with maturity T. Similarly, for each T E T and 
D E V let IV D,T and IVd,t{Q) denote the observed and model implied volatilities with the given 
D,T. The objective function to be minimized is a sum of squared relative errors. We introduce a 
weight factor between the CDS and IV terms to offset a natural overweighting stemming from the 
large number of equity securities relative to the credit securities. Without this factor, the IV terms 
would dominate the calibration and wash out the credit effects we aim to capture. Thus we define the 
objective function to be 

..p.. _ \CDST-CDSTm'' , ^ V- \IVd,t-IVdAQ)? 

'^V^J 2-^ . .2 + (J2 2-^ - — 2 

Tef CDS J' T£T,D£V ^^D,T 

where the subjective value = 7 is chosen to provide a nice balance between the credit and equity 
datasets. The model calibration is required to minimize J over the domain 6 G M x M^: 

6 = argmin0jl(6) 

Remark 3. The above objective function corresponds to the measurement hypothesis that individual 
spreads and implied volatilities are observed with independent Gaussian relative errors. Of course, 
this is likely far from true since in reality, the raw data has been "cleaned up" and transformed by 
Bloomberg in many different ways. Another important point is that the weight = 7 in our objective 
function was chosen by us subjectively to give a nice balance between the credit and equity datasets. 
We tried other two choices of weighting schemes, = 1 and a version weighted by squared bid/ask 
spreads. Both these schemes turn out to under- weight CDS data resulting in a good fit to equity 
options while ruining the fit to CDS spreads. 

The calibration was implemented on a laptop using standard MATLAB. Minimization of J was 
implemented using a quasi-Newton optimization algorithm within MATLAB's function "fmincon": 
typically around 200 evaluations of J were needed to find an acceptable minimum. The model CDS 
and call option prices entering J are computed by one and two dimensional Fast Fourier Transforms 
using MATLAB's functions "FFT" and "FFT2". Since each FFT2 is computationally intensive, and 
calibration involves a great number of evaluations, it was effective to use interpolation from a single 
grid of FFT-generated prices to compute the range of option prices across moneyness. We used 
the truncation and discretization error analysis in [|T3l to optimize these computations. During the 



calibration, we found that the down-and-in barrier option terms in ([29]) are always much smaller than 
the vanilla terms. This can be explained because the credit quality of Ford is reasonably high, and 
therefore the linear transformation {vt, dt)' = R{vt, dt)' generates a point Vt < dt equivalent to a deep 
out-of-the-money vanilla spread option. 



5.3 Calibration Results 

The calibration results for the three models and two dates are shown in Table[T] and Figures 1-9 show 
the model and observed CDS and implied volatility curves in each of these cases. We see that both 
the EXP and VG models lead to CDS and IV curves that capture the correct levels and quantitative 
features of the market data. On the other hand, the GBM hybrid model is able to provide a qualitative, 
not quantitative, fit to the observed market data on the early date, but fails to capture the curves 
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observed on the later date. As a stochastic volatility model, and unlike the Black-Scholes model, it 
does generate an implied volatility smile, but the details of the fitted shapes are not very good. 

One important point is that both LSBM models (but not the GBM model) on both dates lead to 
the implied recovery value i? = 0, at the boundary of the parameter domain and counter to common 
sense. Since reducing R to zero raises the level of CDS curves while leaving IV curves unchanged, 
this observation suggests a mismatch of risk neutral pricing between equity and credit markets. If 
the observed CDS spreads are somewhat high relative to the level of the implied volatility surface, 
the calibration will be forced to choose a zero recovery rate in order to narrow the gap as much as 
possible. In our calibrations it appears that this gap is completely closed with R = 0. A similar 
observation has been made and discussed in Carr and Wu who find a zero recovery rate for 6 out 
of 8 companies studied, leaving a significant gap between model prices and market quotes. 

Table |2] records the balance sheet entries V, D, S implied by each of these calibrations and, for 
comparison, the summary asset and liability values taken from Ford's quarterly reports on the dates 
closest to the calibration dates. Clearly there is no obvious direct correspondence between the balance 
sheet and market values (which appear to be quite sensitive to the choice of model). It will be a 
nontrivial task for the future to investigate the extent to which the balance sheet numbers can be 
predicted by the market implied values. 

6 Conclusions 

We have shown that simple capital structure models that treat the firm value and debt value as cor- 
related stochastic processes can simultaneously capture key empirical features of prices observed in 
credit and equity markets. This framework is a significant extension of the standard structural credit 
risk framework that treats the firm's debt as deterministic or constant over time. We have shown 
that the necessary pricing computations for vanilla equity options and CDS spreads reduce to explicit 
one and two dimensional Fourier transforms, and so implementation of such models is comparable 
in difficulty to Levy process models now used routinely by banks. These promising results indicate 
that it will be worthwhile exploring the performance these models on a variety of firms from different 
countries and sectors. It will also be of interest to explore improvements that may come by using 
other types of time change. It is in principle quite straightforward to move beyond the present TCBM 
model with Levy subordinators to a far broader modeling framework of time-changed Brownian mo- 
tions. A range of implementation and calibration issues will be encountered when considering general 
TCBMs, notably stochastic volatility and jump effects, and addressing such issues will make this a 
major research and development undertaking. 
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2010/07/14 


2011/02/16 






0.2433 


0.2005 






0.1344 


0.1473 




p 


-0.0699 


-0.0143 




b 


0.4966 


0.6948 


VG Model 


c 


0.0474 


0.0240 




R 










V 


3.1796 


3.3393 




d 


2.5036 


2.4973 




RMSE 


0.0804 


0.0271 




- V 


0.2502 


0.2011 






0.1324 


0.1553 




p 


-0.1687 


-0.0383 




b 


0.3700 


0.7232 


Exponential Model 


c 


0.0519 


0.0416 




R 










V 


3.2786 


3.3248 




d 


2.6898 


2.4633 




RMSE 


0.0801 


0.0265 






0.0469 


0.0612 






0.0130 


0.0095 




P 


-0.8175 


-0.9508 


Black-Cox Model 


R 


0.1900 


0.4225 




V 


4.5640 


4.4767 




d 


4.4327 


4.2752 




RMSE 


0.0983 


0.1461 



Table 1: Parameter estimation results and related statistics for the VG, EXP and Black-Cox models. 
"RMSE" denotes the root mean square relative error and is given by (JI(6))^/^ with = 1. 
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FS 


VG 


EXP 


GBM 


Asset 


119.0/109.5 


82.9/107.2 


91.6/105.6 


331.1/334.2 


Debt 


i%3ni.i 


42.2/46.2 


50.8/44.6 


290.3/273.2 



Table 2: Ford's asset and debt values (in $ Bn) reported in the nearest quarterly financial statements 
(June 2010 and March 2011) and estimated from models on July 14 2010 and February 16 2011. 
On those two dates, the number of Ford shares outstanding was approximately 3450 MM shares and 
3800 MM shares respectively. From the financial statements, we take the asset value V to be the total 
current assets plus half of the total long-term assets, and the debt value D to be the current liabilities. 
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Figure 1: CDS market data ("x") versus the GBM model data ("o") on July 14 2010 (left) and 
February 16 20 11 (right). 
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Figure 2: Implied volatility market data ("x") versus the GBM model data ("o") on July 14 2010. 
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Figure 3: Implied volatility market data ("x") versus the GBM model data("o") on February 16 201 1. 
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Figure 4: CDS market data ("x") versus the VG model data ("o") on July 14 2010 (left) and February 
16 2011(right). 
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Figure 5: Implied volatility market data ("x") versus the VG model data ("o") on July 14 2010. 
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Figure 6: Implied volatility market data ("x") versus the VG model data ("o") on February 16 2011. 
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Figure?: CDS market data ("x") versus the EXP model data ("o") on July 14 2010(left) and February 
16 2011(right). 
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Figure 8: Implied volatility market data ("x") versus the EXP model data ("o") on July 14 2010. 
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Figure 9: Implied volatility market data ("x") versus the EXP model data ("o") on February 16 201 1. 
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